Spin Liquid Ground State in the Frustrated Kagome Antiferromagnet MgCu 3 (OH) 6 Cl2 
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We report /iSR experiments on Mg I Cu4_ a! (OH)6Cl2 with x ~ 1, a new material isostructural to 
Herbertsmithite exhibiting regular kagome planes of spin | (Cu 2+ ), and therefore a candidate for 
a spin liquid ground state. We evidence the absence of any magnetic ordering down to 20 mK (~ 
J/10 4 ). We investigate in detail the spin dynamics on well characterized samples in zero and applied 
longitudinal fields and propose a low T defect based interpretation to explain the unconventional 
dynamics observed in the quantum spin liquid phase. 

PACS numbers: 75.10.Kt, 76.75.+i, 75.30.Hx 



Quantum magnetism in frustrated networks has been 
acknowledged for a long time as the ideal playground 
to stabilize new quantum phases. Since the proposal of 
an RVB state by P.W. Anderson [![, advances in the- 
ory has led to the emergence of a rich variety of spin 
liquidphases, including algebraic and gapped spin liq- 
uids [JSj- On the experimental side, Herbertsmithite, 
ZnCu3(OH)gCl2, is among the best materials to explore 
such quantum states, as it combines the highly frustrated 
two-dimensional kagome lattice and quantum spins i of 
Cu 2+ Q. No magnetic freezing has been detected down 
to 50 mK - J/4000 d i] and spin-spin correlations are 
found to be short ranged and dynamical as expected for 
a liquid state Q ■ 

Nevertheless, this material deviates from the "perfect" 
image since it has sizable level of Cu-Zn intersitc mixing. 
The replacement of non-magnetic Zn 2+ ions at the inter- 
layer position by magnetic Cu 2+ causes almost free spin 
| defects that correspond to 5 — 10 % of the total Cu con- 
tent, as determined from various techniques (see Q for a 
review). The exact amount of the complementary defect, 
where a Zn 2+ non-magnetic impurity induces a spin va- 
cancy in the kagome plane, remains difficult to evaluate 
quantitatively, due to similar x-ray scattering factors of 
Cu and Zn and the absence of a straightforward mag- 
netic response. As a result, it varies with samples and 
experimental techniques from a quasi absence of 1 % Q 
(from anomalous x-ray scattering) to a non negligible 5 % 
value [13] (from NMR). 

Very recently, the new series, the Mg-paratacamites 
Mg x Cu 4 _ K (OH) 6 Cl2 isostructural to Zn-paracatamitc, 
has been successfully synthesized. The similar ionic radii 
of Mg 2+ and diamagnctic Zn 2+ ions leads to a min- 
imal difference of the crystal structure and hence ex- 
change pathway, resulting in comparable coupling J ~ 
190 K_[ll|, evaluated by high-temperature series expan- 
sion [la]. If Cu/Mg mixing is still expected, the differ- 
ence in their x-ray scattering factors now enables reli- 



able x-ray structure analysis. Authors of Ref. [l3| suc- 
ceeded to synthesize materials for x < 0.75 where they 
found a ferromagnetic component in macroscopic suscep- 
tibility under Tq = 4 K, attributed to a 3D coupling 
via interlayer Cu 2+ , similar to the Zn case. A minimal 
Mg 2+ substitution within the kagome planes (< 3 %) was 
determined through x-ray diffraction. A different syn- 
thesis was recentlyreported which led to samples with 



0.93 < x < 0.98 [U|, these correspond to the Herbert- 
smithite analogues and correspondingly to model quan- 
tum kagome antiferromagnets. A small ferromagnetic 
fraction at T c ~ 4 — 5 K was tentatively ascribed to an 
impurity phase. 

In this Letter, we first report the /iSR local probe in- 
vestigation down to 20 mK in the MgzCu^^OH^C^ 
Herbertsmithitcs with x <~ 1 which evidences a spin liq- 
uid ground state with no sign of spin freezing. We also 
study the sub-Kelvin spin dynamics of both atacamitcs 
(Mg, Zn) where an x-dependent plateau of \/T\ appears 
with unconventional dynamics. We argue that this re- 
laxation is driven by interlayer Cu 2+ ions and discuss its 
origin. 

The experiments were carried out on Mg atacamite 
powder samples close to the idealized Herbertsmithite 
structure, x = 0.84, 0.92, and 1.21, synthesized following 
a hydrothermal route described elsewhe re , as well as 
on the x = 1 Zn atacamite sample from [l4| . The ratio of 
Cu/Mg and Cu/Zn was determined by ICP-AES (Table 
[I]). Refinements of x-ray diffraction data for the Mg com- 
pounds give separately the Cu (n) and Mg (p) occupan- 
cies of the interlayer and the kagome sites correspond- 
ing to the formula (Cui_ p Mp) 3 (Mi_„Cu„)(OH) 6 Cl2 
(M=Mg, Zn) and thus an independent determination of 
x = 3p— n+1 (Table[T|. yuSR experiments were performed 
at the ISIS and PSI facilities in zero and longitudinal ap- 
plied field configurations down to 20 mK. We also mea- 
sured DC magnetic susceptibility on a standard Quantum 
Design SQUID magnetometer in the 1.8 — 300 K T-range 
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TABLE I: Chemical composition determined through 
ICP, x-ray refinements and saturated magnetization for 
(Cui_ p Mp)3(Mi_ n Cu„)(OH) 6 Cl2 where M=Mg, Zn. Each 
site occupancy is constrained to unity. 



Element M 


Mg 


Mg 


Zn 


Mg 


ICP analysis 


x = 0.84 


x = 0.92 


x = 1 


x = 1.25 


x-ray analysis 


x = 0.83 


x = 0.91 




x = 1.21 


P 


0.04 


0.06 




0.12 


n 


0.29 


0.27 




0.15 


Magnetization 










n 




0.266(4) 


0.217(5) 


0.186(4) 




Field (T) 



FIG. 1: Interlayer Cu 2+ contribution Mdef extracted from 
magnetization measurements and normalized by Af sat = 
Naq^bS — 6143 emu. The saturation value gives n, the 
amount of interlayer Cu 2+ per formula unit. The lines repro- 
duce a Brillouin fit with an interacting scale of 6 = 0.8 K. 
Inset: Normalized magnetization at T = 1.75 K up to 14 T. 
Lines are fits for the linear component. 



and magnetization at 1.75 K up to 14 T with a Cryogenic 
Vibrating Sample Magnetometer. 

In order to reveal the magnetic response of the inter- 
plane defects we measured the magnetization of the Mg 
(x = 0.92 and x = 1.21) and Zn (x = 1) samples at low 
T = 1.75 K and up to 14 T. Following the analysis in 
Ref. [HI i we divide the total magnetization into two con- 
tributions, a Brillouin-likc component arising from the 
quasi-free interlayer Cu 2+ and a linear term coming from 
the strongly coupled Cu 2+ of the kagome planes (Fig[T|) . 
By subtraction of the latter one, we have access to the 
saturated magnetization Md c f = uNa QUb S of the inter- 



layer Cu 2+ . Taking g = 2.2 from Ref. [l5|, we determine 
n which agrees reasonably well with the refinements for 
the Mg samples and provides a value for the Zn sam- 
ple otherwise inaccessible by structural studies (Table [TJ) . 
These results confirm quantitatively the existence and 
the quasi-free behavior of the interplanc defects in the 



Mg and Zn Hcrbertsmithitcs. 

We now turn to the local /iSR investigations. /iSR 
is very sensitive to any small magnetic field (down to 
~ 0.1 G) and is therefore a powerful tool to detect any 
frozen moment. The muon is implanted inside the vol- 
ume and, as a positive charged particle, will stop in the 
vicinity of a negative environment, i.e. either near OH~ 
(80 %) or CI" (20 %). The polarizations from (ZF) 
to 2500 G longitudinal applied field (LF) are reported in 
Fig|2] Our experiments demonstrate a similar magnetic 
behavior of the Mg and Zn compounds, as expected due 
to their crystallographic similarities. Following the for- 
mer work of Ref. the ZF polarization is fitted on a 
high statistics run at 50 mK by P(t) = F nuc i(t)e _(At) ' 9 . 
-Pnuci(£) depends only on static fields from surrounding 
nuclei and A stands for a small dynamical relaxation. The 
well defined oscillation of P n uci(i) is due to the formation 
of a /i-O-H complex [l6| . and from the field experienced 
by the muon, if M _oH = 7.8 G, one estimates a distance 



of 1.5 A to hydrogen 17 1. The ZF polarization is found 
identical in the extended T-rangc 0.05 — 20 K, except for 
a slight variation of A. Therefore, wc conclude that there 
is no magnetic ordering of the electronic spins down to 
T = 20 mK (~ J/10 4 ) for x ~ 1, as in the Zn counter- 
part. 

The spin dynamics are revealed by experiments under 
longitudinal fields. From 10 to 100 G (> 10# m _ O h), the 
static relaxation of nuclear origin is progressively decou- 
pled as expected (FiglJ). Above 100 G, P n uci(*) = 1 and 
a single stretched exponential fit of P(t) yields the re- 
laxation rate A = 1/T-f of electronic origin. A = \jTy 
is linked to the spin-spin time correlation function by 
1 / T i ~ Jo +00 ( S ( t ) S (°)) cos( lfl H LF t)dt, H LF is the longi- 
tudinal applied magnetic field. Starting from a high-T 
value ~ 0.005 /is -1 , A increases upon cooling below 1 K 
and saturates at lower temperature at a value which dif- 
fers from sample to sample (FigJ3]). This increase of A 
indicates a dynamical slowing down. We notice that the 
shape of the relaxation at low T slightly changes between 
Mg (/3 = 0.9) and Zn compounds ((3 = 1.1). The origin 
could be related to a presently unclear difference between 
the spin-spin time correlation in the two materials. 

We now argue that the muon is dominantly coupled 
to the interlayer Cu 2+ moments from the following ex- 
perimental observations: (i) The low-T muon shift K* 1 
was shown to track the intersite defect susceptibility [181 . 
This is in fair agreement with recent 2 D NMR study [19( 
but contrasts with NMR results on 17 which is strongly 
coupled to the planes, (ii) From Fig[3] (b) and the in- 
spection of Table [l] A is found to increase linearly with 
n, the concentration of intersite defects, but is anticor- 
related with the level of in plane defects, p. Since the 
distance between Zn sites is d = 6.12 A, the linearity in 
n can be easily explained by an "all or nothing" model: 
some muons (n) sit next to an interlayer Cu 2+ defect 
with a relaxation rate Ai whereas the others (1 — n) stay 
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FIG. 2: Upper panel: Polarization in zero field (ZF, black 
circles) and under longitudinal applied fields (LF) from 10 to 
2500 G. The black line is a fit (see text). Bottom panel: The 
relaxation is faster at low temperatures and when n increases. 
Lines are stretched exponential fits. 



FIG. 3: (a) T-dependence of relaxation rates A for different 
samples, (b) Plot of the plateau value of A = 1/21 versus n. 
The black line is a linear fit. The error bars on n account for 
the dispersion between M(H) and x-ray characterization. 



far from a defect with a relaxation rate A2 — >• 0. Under 
the valid condition A1.2 t<l, we get a total polarization 
P(t) = ncxp(— Ait) + (1 — n)exp(— X 2 t) ~ exp(— nXit), 
which explains the linear n-dependence of A. (Hi) From 
the scaling of versus ^ bulk 5 one can extract the cou- 
pling constant = 0.08 T//_tg (l8| . consistent with a 
dipolar interaction. From 17 NMR Tf 7 , the contri- 
bution of the kagome planes Cu 2+ to the muon relax- 
ation can then be estimated using: l/Tf = l/Tf 7 x 
(7^/7," A^) 2 - 3 • 10~ 5 [is- 1 . This leads at high-T to 
a value 150 times smaller than the measured one, which 
calls for a different source of relaxation, (iv) Finally in 
the high-T Moriya paramagnetic limit , 1 /Tf is given 
by l/Tf = 2 7 2 iT 2 n/^ with v = ^/AJ l2 zS{S + l)/3nh 2 
and Hfj, = gA^S(S+l)/3 



constant value l/Tf ~ 5 • 10~ 3 



880 G. From the high-T 
/its , an average value 
?i ~ 0.22 and a number of nearest neighbors z = 6, one 
obtains the coupling between an interlayer Cu 2+ and a 
kagome Cu 2+ J' ~ 3 K, in rough agreement with the 
temperature of the magnetic ordering Tq = 6 K in the 



clinoatacamite Cu2(OH)3Cl parent compound [211.122 1. 
In this context, the slowing down below 1 K of intcrsitc 



defects can be attributed either to the strengthening of 
the correlations with the two nearby kagome planes or to 
a coupling between two intersite defects mediated by the 
kagome plane to which they are coupled. The tempera- 
ture scale of 1 K is in line with subtle effects detected in 
magnetization measurements flij . 

We now discuss the T — > value of the relaxation rate. 
In order to reveal the dynamics of the correlated regime 
under 1 K, we probe the excitation spectrum S(uj) in the 
low energy range by applying a longitudinal field TJlf = 
w/7^, TJlf < 0.25 T. Two scenarios can be considered 
based on different correlation functions S(t) = (S(/j)S(0)) 
for the interplane spins. 

Exponential correlation function S(t) = e~ yt associated 
with field induced polarization of the interlayer defects. 
Such correlation leads to the usual lorentzian spectral 
density S(u>) and 



Ai = 



v 2 + llH- 



(1) 



LF 



where -f/fl uc t is the fluctuating component of the field at 
the muon site perpendicular to its initial polarization, v 
is the fluctuation frequency, and 7 M = 851.6 Mrad/s/T 
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is the muon gyromagnetic factor. However, not only the 
expected i/ LF dependence of A is not convincing (Fig|4}, 
but a forced fit will result in an unphysical Hfi uct ~ 20 G, 
while the lowest possible dipolar field value is 200 G 
which corresponds to a maximal distance of the fi + to 
intersite Cu 2+ . We therefore propose that the fluctu- 
ating field is reduced when the S = interlayer starts 
to be polarized in the external applied field H^p. In a 
mean field Brillouin approach, the fluctuating moment is 
m fluct = MB (i_tanh(. 9 ^ B Si/ LF /fc B (T+(9))) where is in- 
troduced to account for interactions. The reduced value 
of -Hfluct = rrfi^H^/fiB can be injected in EqQ] and 
leads to only two shared fit parameters to account for the 
T^(Hlf) variation for all x. Although the obtained fits 
are not perfect, this approach yields reasonable physical 
parameters v = 100±20 GHz and 9 = 0.7±0.2 K that are 
consistent with both magnetization measurements and 
the 1 K T-scale where the slowing down occurs. 
Power law correlation function S(t) = (l/t) 1- " 
The corresponding spectral density S(u>) yields a power 
law relation T-f oc Lo a . This more exotic approach is at 
play (a = 0.35) in the spin liquid case of the S = ^ an- 
tifcrromagnctic chain [23J, where a spinon continuum of 
excitations is well established. Such a spectral density 
(a = 1) was also invoked for the spin liquid pyrochlorc 
TbTi 2 7 dHH]. A perfect fit of our data can be found 
with 1/Ai = T^(x) + A • cj - 63 , which would point to an 
exotic relaxation channel for the intersite defect. Using 
the scaling of V'T Q - (T/w) Q tanh(w//3T) from neutron 
experiments |26ll. this leads to 1/T\ ~ ksTx" '/w ~ u>~ a 
(a = 0.66 [26() by means of the fluctuation dissipation 
theorem, in agreement with the value reported here. This 
points to a common origin and following [27| one could 
invoke a distribution of couplings between intersite de- 
fects as a source of the power law. Another possibility 
could be that the kagome planes dynamics drives that 
of the intersite defects. The inconsistency between the 
T-plateau of the relaxation rate measured at low fields 
and the T~ - 7 dependence of T^ 7 at 7 T would then 
require that the field impacts on the kagome plane dy- 
namics. 

In conclusion, the spin liquid behavior of the new 
kagome antifcrromagnet Mg^Ci^-^OH^C^ for x > 
0.84 is clearly established by our /iSR experiments. Our 
data show that the measured interlayer site dynamics dif- 
fers from those of the kagome plane one at 7 T. This calls 
for a careful inspection of probes which integrate both re- 
sponses and/or the field impact on relaxation. Whether 
the exotic character of the /x + relaxation might also relate 
to the spin liquid behavior of the kagome lattice is still a 
matter of speculation. The ability to refine the structure 
for Mg-Herbertsmithite in a reliable manner opens the 
possibility to control the level of defects and to discrim- 
inate between the various sources of dynamics at low T. 
Our results open routes for future investigations of the 
kagome spin liquid ground state in well controlled mate- 



1 1 1 
" T = 50mK 


1 ' 1 ' 1 ' 1 


* M 9i2, 




■ Zn 1 


















Lorentzian S(e») 




— - Lorentzian S(to) and field effect - 

S{to) ~ CO 



0.10 0.15 0.20 0.25 

H |F (T) 



FIG. 4: The T\ dependence with external field #lf underlines 
the unconventional spin dynamics at T = 50 mK. The lines 
refer to different models (see text). 
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